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ABSTRACT 



Context, Cyclotron resonant scattering features are an essential tool for studying the magnetic field of neutron stars. The fundamental 
line provides a measure of the field strength, while the harmonic lines provide information about the structure and configuration of 
the magnetic field. Until now only a handful of sources are known to display more than one cyclotron line and only two of them have 
shown a series of harmonics. 

Aims. The aim of this work is to see the first harmonic cyclotron line, confirming the fundamental line at ~22 keV, thus increasing the 
number of sources with detected harmonic cyclotron lines. 

Methods. To investigate the presence of absorption or emission lines in the spectra, we have combined RXTE and INTEGRAL spectra. 
We modeled the 3-100 keV continuum emission with a power law with an exponential cut off and look for the second absorption 
feature. 

Results. We show evidence of an unknown cyclotron line at ~47 keV (the first harmonic) in the phase-averaged X-ray spectra of 
4U 1538-52. This line is detected by several telescopes at different epochs, even though the S/N of each individual spectrum is low. 
Conclusions. We conclude that the line-like absorption is a real feature, and the most straightforward interpretation is that it is the 
first harmonic, thus making 4U 1538-52 the fifth X-ray pulsar with more than one cyclotron line. 
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1. Introduction 

Cyclotron resonant scattering features (CRSFs), usually referred 
to as 'cyclotron lines', have proved to be powerful tools for di- 
rectly studying the magnetic field in neutron stars. CRSFs are 
present in the hard X-ray spectra of several X-ray pulsars and 
originate in the "cyclotron process" under extreme conditions. 
Through E cyc = 11.6fi 12 x (1 + z)~ l keV (the '12-B-12 law', 
where z is the gravitational redshift), an energy of the fundamen- 
tal feature in the hard X-rays indicates that the magnetic fields 
are rather strong (B ~ 10 12 G). Under such conditions, the inter- 
action of the electrons and radiation must be treated quantum- 
mechanically. The behaviour of an electron in a strong magnetic 
field implies that the electron energy must be quantized in so- 
called Landau levels. These absorption features stem from the 
resonant scattering of photons by electrons, also referred to as 
cyclotron lines. 

While the fundamental energy of the cyclotron line pro- 
vides valuable information about the magnitude of the field, 
it is only through the detection and the analysis of the har- 
monic lines that we can get direct information about the ge- 
ometrical configuration of the B field (Harding & Daugherty 



1991; |Araya-G6chez & Harding] 2000; ISchonherr etaLl 2007). 
However, to date, only in a handful of systems have harmonic 
lines been discovered, and only two systems have shown more 
than two ( |Santangelo et al.| 1999 & ICoburn et all 2005). It is 
therefore paramount to add as many systems to this selected 
group as we can. 

In this work, we present a spectral analysis of the high mass 
X-ray binary pulsar 4U 1538-52. It is an eclipsing system con- 
sisting of the B0 I supergiant star QV Nor and a neutron star 
with an orbital period of ~3.728 days (Clark 2000). The orbital 
eccentricity is ~0.08 (Corbet et al. 1993 ), although more recently 
a higher value of ~0.17 was deduced by Clark (2000). The X-ray 
eclipse lasts ~0.6 days (Becker et al. |1977t . The system is fairly 
bright in X-rays. The estimated flux is ~ (5 - 20) x 10~ 10 erg 
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s cm ~- in the 3-100keV range (Rodes l2007l) . Thus, assuming 
a distance of the source of ~5.5 kpc (Becker et al. 1 19771 Parkes 
et al. 1 19781 1 and an isotropic emission, the luminosity follows 
~ (2 - 7) x 10 36 erg/s. The magnetized neutron star has a spin 
period of -529 s (Davison [T977l Becker et al. [T977l >. 

The pulse-phase averaged X-ray spectrum of 4U 1538-52 
has usually been described either by an absorbed power law 
modified by a high energy cutoff, a power law modified by a 
Fermi-Dirac cutoff, or by two power laws with indices of oppo- 
site sign multiplied by an exponential cutoff (the NPEX model, 
Mihara 1995; Rodes et al . 120061 . In addition to these continuum 
models, an iron fluorescence line at ~6.4 keV and a cyclotron 
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Fig. 1. Combined spectrum and model of data obtained with PCA (3-20 keV) and HEXTE (17-100 keV). Both data sets belong to 
the run carried out in 2001 and their orbital phases are 0.53 and 0.66, respectively. Bottom panels show the residuals in units of <x 
with respect to the model (see Section lXTl for details). 



resonant scattering feature at ~20 keV discovered by Ginga 
(Clark et al. 1990 ) are needed to describe the data. The variabil- 
ity of this CRSF was studied by Rodes-Roca et al. ( 120081 1. Rossi 
X-ray Timing Explorer (RXTE) (Coburn [200Tb and BeppoSAX 
data (Robba et al. 2001) did not show evidence of the first har- 
monic at ~40 keV. Robba et al. ( 12001 1 ) presented some evidence 
of an absorption feature around 50 keV; however, because of the 
lack of a signal-to-noise ratio of the spectrum at these energies, 
the feature could not be confirmed. 

In this paper, we report on the 3-100 keV analysis based on 
the observations of 4U 1538-52 performed by the RXTE and 
INTEGRAL satellites. In Sect. [2] we describe the observations 
and data analysis. In Sect. [3] spectral analysis are presented, and 
summarized in Sect. [4] 

2. Observations 

2.1. RXTE data 

To study the presence of spectral features, we have used all 
archival data from RXTE on this source: three observations have 
been carried out in 1996, 1997, and 2001. The first one, is 
a monthly observation between 1996 November 24 and 1997 
December 13. The second and third ones cover a complete or- 
bital period. In our analysis we used data from both RXTE point- 
ing instruments, the Proportional Counter Array (PCA) and the 
High Energy X-ray Timing Experiment (HEXTE). 

To extract the spectra, we used the standard RXTE analysis 
software FTOOLSQ This package takes care of the modeling 
of PCA background, the dead time corrections of HEXTE data 
and generates the appropriate response matrices for the spectral 
analysis. 

The PCA consists of five co-aligned Xenon proportional 
counter units with a total effective area of -6000 cm 2 and a 
nominal energy range from 2 keV to over 60 keV (Jahoda et 
al. [1996 ). However, because of response problems above ~20 
keV and the Xenon-K edge around 30 keV, we restricted the use 
of the PCA to the energy range from 3 keV to 20 keV (see also 
Kreykenbohm et al. 120021 ). Systematic uncertainties are taken 
into account by the standard spectral analysis. 
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The HEXTE consists of two clusters of four NaI(Tl)/CsI(Na) 
Phoswich scintillation detectors with a total net detector area 
of 1600 cm 2 . These detectors are sensitive from 15 keV to 
250 keV (Rotschild et al. 119981 ), however, response matrix, in- 
strument background and source count rate, limit the energy 
range from 17 to 100 keV. Background subtraction in HEXTE is 
done by source-background swapping of the two clusters every 
32 s throughout the observation. For the HEXTE, the response 
matrices were generated with HXTRSP, version 3.1. We used 
HXTDEAD version 2.0.0 to correct for the dead time. In order 
to improve the statistical significance of the data, we added the 
data of both HEXTE clusters and created an appropriate response 
matrix by using a 1:0.75 weighting to account for the loss of a 
detector in the second cluster. We also binned several channels 
together of the HEXTE data at higher energies and chose the bin- 
ning as a compromise between increased statistical significance 
while retaining a reasonable energy resolution. 

2.2. INTEGRAL data 

INTEGRAL [(Winkler et al. 2003)| has a unique broad band capa- 
bility thanks to its four science instruments (imager IBIS, spec- 
trometer SPI, X-ray monitor JEM-X, and optical monitor OMC) 
which allow us to study a source from 3 keV up to 10 MeV and in 
the optical simultaneously. The imager IBIS has a very large field 
of view of 9° x 9° which allows us to observe many sources at 
the same time. Together with a large collecting area of 2600 cm 2 
and decent energy resolution of 9% at 100 keV makes IBIS the 
prime instrument for our analysis. 

Since its launch on October 17 2002, INTEGRAL observa- 
tory has been constantly collecting a wealth of data. INTEGRAL 
data are organized in revolutions (i.e. 72 h long satellite orbits 
around the Earth) and then science windows which are typically 
1800 s to 3600 s long. 

The source was observed by the hard X-ray imager (IBIS) 
camera on board INTEGRAL during the regular scans of 
the Galactic plane and the Norma survey. Fig. [2] shows an 
IBIS/ISGRI mosaic image of the Norma Arm region in the en- 
ergy range 20-100 keV. The source 4U 1538-52 (labelled as 
H 1538-522) is clearly detected and its position is well deter- 
mined, allowing us to extract its spectrum without contamination 
of other sources in the field of view of the instrument. 
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Fig. 2. IBIS/ISGRI mosaic image of the Norma Arm region 
where the source is located and clearly detected. 



Table 1. Details of observations. 
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Fig. 3. Combined spectrum and model obtained with JEM-X 
(left), ISGRI (upper right), and SPI (lower right). The contin- 
uum is modeled by the cutoffpl model without any cyclotron 
lines applied. The bottom panel shows the residuals in units of 
cr with respect to the model. The first harmonic at ~47 keV is 
already evident in the raw data. 



Satellite/Instrument 


MJD 


Exposure 


Orbital 






(ks) 


phase 


RXTE/PCA 


51925.25 


9.296 


0.53 


RXTE/HEXTE 


51925.25 


6.418 


0.53 


RXTE/PCA 


51925.74 


12.592 


0.66 


RXTE/HEXTE 


51925.74 


8.261 


0.66 


INTEGRAL/JEM-X 1 


53407.7-53826.3 


360 




INTEGRAL/ISGRI 


52651.4-53756.2 


720 




INTEGRAL/SPI 


52651.4-53941.4 


1907 





To extract INTEGRAL data, we used the Offline Science 
Analysis Software (OSA) version 7.0, as distributed by the 
Integral Science Data Centre (ISDC) following the respective 
cookbook instructions. For IBIS, we selected all public data up 
to revolution 400 with the source in the fully coded field of view 
of the instrument, resulting in ~400 Science Windows (ScWs). 
In a first step we created a mosaic of the full data set to ob- 
tain a catalog of the detected sources in the field of view. Since 
4U 1538-522 is relatively close to the galactic centre, the cata- 
log contains other 15 sources, among them several bright sources 
such as 4U 1700-377. We then used this catalog to extract the 
spectrum of 4U 1538-522 using all 400 ScWs in order to ob- 
tain a high signal-to-noise ratio. Since all other sources in the 
field are several degrees away, source confusion is not a prob- 
lem. For JEM-X we used all data of monitor number one (JEM-X 
1) within 2° to ensure a reliable spectrum resulting in approxi- 
mately 55 ScWs. And for SPI, the selection of all available data 
within a 8° radius (the fully coded field of view of the instru- 
ment) allowed us to obtain 969 ScWs. The selected spectra then 
have an effective exposure time of 720ksec for IBIS, 1907 ksec 
for SPI, and 360 ksec for JEM-X. 

See Table Q] for the details of all observations which have 
been used in this work. 



3. Spectral analysis 

For spectral analysis we used the XSPEC dArnaud 19 96) fitting 
package, released as a part of XANADU in the HEASoft tools. 

3.1. RXTE analysis 

Fig. Q] shows PCA and HEXTE phase-averaged spectrum at or- 
bital phases 0.53 (left panel) and 0.66 (right panel). The raw 
data together the best-fit model and the residuals of the fit as 
the difference between observed flux and model flux divided by 
the uncertainty of the observed flux, i. e. in units of cr, are in- 
cluded in this plot. The dip of the second cyclotron line at ~ 
47 keV is apparent in the raw data. The RXTE continuum is 
properly described by an absorbed powerlaw modified by an 
exponential cutoff at E cut = 9.4+iL keV (cutoffpl in XSPEC 
dArnaud 19 96)). see Table [2] Analytically, cutoffpl is given by 
the equation, 

cutoff pl(E) = K (E/lkeVy a <T £ / £ »' , (1) 

where a is the power law photon index and E cut is the cutoff en- 
ergy in keV. An iron emission line at 6.4 keV was also added. 
The proper description of the continuum, however, further re- 
quired a very broad Gaussian emission line component at 12 keV 
to account for a flux excess at these energies. This phenomeno- 
logical component has been used to describe the continuum of 
several X-ray pulsars (Coburn et al. [2002 concluded that this 
bump is an inherent feature in the spectra of accreting pulsars). 
Although it has no clear physical interpretation, its use here is 
justified as we concentrate on the presence of absorption fea- 
tures at significantly higher energies. See Table[2]for the best fit 
spectral parameters for the data set in Fig.Q] Data from PCA and 
HEXTE have been fitted simultaneously. All uncertainties refer 
to a single parameter at the 90% (Ay 2 = 2.71) confidence limifl 
A significant absorption feature is present in the residuals 
around ~47 keV. When other models are used to describe the 
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Table 2. Fitted parameters for the RXTE spectra in Fig.Q] 



Table 4. Fitted parameters for the fundamental and the first har- 
monic CRSF. 



Component 


Parameter 


Left panel 


Right panel 


Continuum 


a 

E cut (keV) 
N H (10 22 cm" 2 ) 


-0,14 

11.4+2° 
1 +0 6 


0.52±0.11 
9.4 + 'i 

-0.23 

1 05 + 010 

±." -0.20 


Gaussian 
emission line 


Eh,,,,,,, (keV) 
o-bump (keV) 


12 54 +0 2 ' 

• J -0=1 
o 19+0.22 

-0.23 


12 15 +01 ° 

-0.Q9 
9 Q9+0.08 


Fluorescence 
iron line 


E FeK (keV) 
cr FeK (keV) 


6 54 +014 

"•^-0.19 


6 56 +0 03 
0.1 (frozen) 






1.2(51) 


4.4(52) 



Parameter 


SPI 


ISGRI 


HEXTE 


Fundamental 








<t ci (keV) 
1" Harmonic 


01 '7+0.9 
ZZ -'_2 4 

0.8*g 

a +4 
J -1.7 


99 n+i.o 
5 +1 6 

-0.18 

3 +5 


9f\ 7+1.3 

5 +1 - 3 

— 0.13 

4+6 


E C1 (keV) 
a„ (keV) 


42 +4 
9+1.0 

Q + ll 

-1.9 


47+3 
0.8+ 021 


48+' 
OKI 


^(dof) 


1.0(14) 


1.1(18) 


0.8(16) 



Table 3. Fitted parameters for the INTEGRAL spectrum in Fig. [3] 



Component 


Parameter 




Continuum 


a 


0.29+0.07 




E m (keV) 


9-061 22 


Normalization 


CjEM-X 


0.90±0.03 




C ISGRI 


1.0 (frozen) 




Cspi 


i 9Q +0.08 
1,z - 7 -0.07 






5.5(68) 



underlying continuum, i.e. an absorbed powerlaw with a Fermi- 
Dirac cutoff dTanakal 1986) or the NPEX model (Mihara |l995V 
the overall residuals were higher, but none of them was able to 
account for the absorption feature at ~47 keV. 

3.2. INTEGRAL analysis 

In Fig. [3] we present the combined JEM-X, ISGRI and SPI spec- 
trum using all data. The spectra of all three instruments were 
fitted simultaneously with a cutoffpl model, using the latest re- 
sponse matrices provided by the OS A software (see Section 12721 
for more details). A factor term was included in the model to 
allow for the adjustement of efficiencies between different in- 
struments. No photoelectric absorption has been used since the 
INTEGRAL spectra start only at 4 ke V and none was required by 
the data. An absorption column below 10 23 cm -2 is expected to 
have no noticeable effect below 5 keV. Nevertheless JEM-X data 
are needed to constrain the X-ray continuum model. Likewise, 
neither the very broad Gaussian at 12 keV, nor the iron fluores- 
cence line at 6.4 keV are required by the JEM-X data and were 
therefore not added to the model. As can clearly be seen in Fig. [3] 
significant absorption line structures are present at ~20keV (the 
well known CRSF) and ~ 47keV after fitting the continuum 
with a cutoffpl model (best-fit parameters of continuum are given 
in Table [3). Data from JEM-X, ISGRI and SPI have been fitted 
simultaneously. 

3.3. Combined RXTE and INTEGRAL analysis 

To further explore this second feature and in order to achieve 
the highest significance at high energies, we jointly fitted the 
data of the high energy instruments of the two satellites, HEXTE, 
ISGRI, and SPI. We fixed the continuum parameters a, the power 
law photon index and E cut , the cutoff energy of the exponential 
cutoff, to the values given in Table [3] 



In Table|4]we show the best fit parameters for each spectrum 
we used in Fig. [4] Data from SPI, ISGRI and HEXTE have been 
fitted individually. As the continuum parameters are frozen, we 
have not included them in this Table. The results of the spec- 
tral fits from three different instruments show that the values de- 
scribing the shape of the cyclotron lines are consistent with one 
to another within uncertainties. 




20 50 
channel energy (keV) 



Fig. 4. Spectra and continuum model (cutoffpl) obtained with 
ISGRI, SPI {lower), and HEXTE. The bottom panel shows the 
residuals in units of cr of a fit without including a cyclotron line. 

As evident from Fig. @] the two cyclotron absorption fea- 
tures are clearly seen in the raw data of the three instruments. 
Although SPI has a very high energy resolution, the compared to 
ISGRI rather low signal-to-noise ratio required the use of rather 
broad energy bins. Therefore we will use ISGRI as the prime in- 
strument in our study together with HEXTE while SPI data will 
be used for comparison. In summary, to get the best fit param- 
eters of the absorption features, we have combined ISGRI and 
HEXTE data including a factor for free normalization between 
the two instruments (see Fig. [5]top panel). The two absorption 
features at ~22 keV and -47 keV are modeled using the cyclabs 
model from XSPEC dArnaud 19961 

We started by modeling the fundamental cyclotron line that 
was discovered by the Ginga satellite observatory (Clark et 
al. 1 19901 1 (see Fig. [5] second panel). After the inclusion of the 
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Table 5. Best fit parameters for the fundamental and the first 
harmonic CRSF. 



Component 


Parameter 






Continuum 


a 


(keV) 


0.29 (frozen) 
9.1 (frozen) 


r UlllLilllL 1 Hill 


17 

Z C , 

T* 


(\ca\r\ 
(,kcv ; 

(keV) 


91 4 +M 

43 +0 3 

"•^ -0.07 
o+4 

3 -2 


I st Harmonic 


E C2 

°"C2 


(keV) 
(keV) 


47 r 2 - 2 

°- 88 -024 

5 +5 

-2.3 

1.4(41) 



fundamental cyclotron line at £ ~ 21.4 keV, the xl improves 
from 5.5 for 47 degrees of freedom (dof) to 3.4, for 44 dof (F- 
test: 3.0 x 10~ 5 ). Although residuals improve significantly (see 
Fig.|5]third panel), another absorption feature can be seen at ~47 
keV. Including a second cyclotron line at ~ 47 keV improves the 
fit further resulting in axl of 1.4 for 41 dof (F-test: 1.7xl0~ 8 see 
Fig Hbottom panel). ISGRI data (filled grey circles in Fig.[5]l has 
a far better resolution than HEXTE at these energies. However, 
the overall shape is virtually identical, ruling out any instrumen- 
tal or circumstantial effects. The final fit parameters are given in 
Table [5] Data from HEXTE and ISGRI have been fitted simulta- 
neously. 

The F-test is known to be problematic when used to test the 
significance of an additional spectral feature (see Protassov & 
van Dik 2002), even if systematic uncertainties are not an issue. 
However, the low false alarm probabilities may make the detec- 
tion of the line stable against even crude mistakes in the compu- 
tation of the significance (Krevkenbohm l2004l l . Therefore, tak- 
ing into account these caveats, we can conclude that the the first 
harmonic CRSF is detected with high significance in the spec- 
trum of this source. 

Although the uncertainties of the determined values are 
rather large, these values are consistent with one another within 
uncertainties. This makes it extremely unlikely that the ~47 keV 
feature results from a calibration problem. Nevertheless, the line 
parameters depend slightly on the shape of the continuum (for 
example, using the NPEX component (Mihara 1995) and two 
cyclotron absorption lines, we obtained the line centre energies 
at ~22 keV and ~44 keV, for the fundamental and first harmonic 
respectively). 

In short, this feature has been found to be present under the 
following circumstances: 

- Three different telescopes and instruments: Narrow Field 
Instruments (NFIs) on board BeppoSAX (Robba et al. 
I200TT) . HEXTE on board RXTE, and IBIS/ISGRI on board 
INTEGRAL (this work). 

- Different epochs within the same telescope (i.e. 1996, 1997, 
2001 using RXTE). 

- Different orbital phases within a given epoch and instrument 
(Fig.0. 

We therefore conclude, that despite the relatively low signifi- 
cance of the feature in an individual observation, its detection in 
different epochs, configuration of the binary system, and differ- 
ent instrumentation, makes this feature highly significant. The 
most direct explanation for this absorption line like feature is 
that it is the first harmonic cyclotron line. 
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Fig. 5. Cyclotron line modeling of the phase-averaged spectra 
HEXTE/ISGRI of 4U 1538-52. Top panel: Spectra and best fit 
model (cutoffpl and two cyclotron lines) obtained with HEXTE 
and ISGRI. Bottom panels show the residuals in units of cr 
for models taking different numbers of cyclotron lines into 
account. Second panel: Without cyclotron lines. Third panel: 
Fundamental cyclotron line at ~22 keV. Bottom panel: Two cy- 
clotron lines, the fundamental at ~22 keV and the first harmonic 
~47 keV (see Table|5]for the best fit values). 

4. Summary and discussion 

We presented the spectral analysis of 4U 1538-52 using data 
from RXTE and INTEGRAL. We present evidence for a previ- 
ously unknown absorption line like feature in the phase-averaged 
spectrum of the source. As we have shown in Section[3] we have 
been able to achieve a good fit to the phase-averaged spectra by 
including a Lorentzian absorption line at ~47 ke V into the model 
(see Fig . O ■ This absorption line is clearly visible whenever the 
signal-to-noise ratio in the spectrum is good enough to allow an 
analysis of the data. The most straightforward interpretation for 
this feature is that it is the first harmonic of the ~ 22 keV funda- 
mental CRSF. 

According to the theory, cyclotron lines are due to the res- 
onant scattering of photons by electrons whose energies are 
quantized into Landau levels by the magnetic field (Meszaros 
119921 ). The quantized energy levels of the electrons are harmon- 
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ically spaced in the first order, such that the first harmonic line 
should be placed at twice the energy of the fundamental line, i.e. 
2 x E cyc x 42 keV. In reality, however, the coupling factor be- 
tween the fundamental and first harmonic is with ~ 2.20 slightly 
higher than 2.0. This anharmonic spacing, however, has been 
observed already in several systems where more than one line is 
present. As explained by Schonherret al. (2007), the relativistic 
photon-electron scattering already produces some anharmonic- 
ity, because photons with energies close to the Landau levels 
may not escape the plasma if their energies are not changed by 
inelastic scattering. This, however, can not be the only reason be- 
cause some systems show an anharmonic spacing larger than that 
predicted by this effect. A possibility to explain this difference 
is to take into account that the optical depths of the fundamental 
and the first harmonic could be different if they are formed at dif- 
ferent heights above the neutron star. With increasing height, the 
strength of the magnetic field decreases resulting in a different 
CRSF energy. Another possibility is to consider a displacement 
of the magnetic dipole which would also explain the difference 
of energy of the two lines if the lines originate from the different 
poles of the neutron star. Therefore a significant phase depen- 
dence of the strength of the both lines is expected, however, the 
low signal-to-noise ratio at higher energies prevents us to test 
this hypothesis with the current data sets. 
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